IMPULSE CONTROL OF MULTIDIMENSIONAL JUMP DIFFUSIONS 



MARK H. A. DAVIS*, XIN GUO+, AND GUOLIANG WU* 

Abstract. This paper studies regularity property of the value function for an infinite-horizon discounted cost impulse 
control problem, where the underlying controlled process is a multidimensional jump diffusion with possibly 'infinite-activity' 
jumps. Surprisingly, despite these jumps, we obtain the same degree of regularity as for the diffusion case, at least when the 
jump satisfies certain intcgrability conditions. 

\ 1. Introduction. This paper is concerned with regularity of the value function in an impulse control 

■ problem for an n-dimensional jump diffusion process X(t). 

In the absence of control, the stochastic process X(t) is governed by the following SDE: 

(N ' 

O; dX(t) = n(X(t-))dt + a(X(t-))dW(t)+ I j{X{r),z)N{dt,dz), X(0) = x. (1.1) 



Here W(t) is an m-dimensional Brownian motion and N(-, •) a Poisson random measure on R + x M 1 , with 
W and N independent. The Levy measure := E(A r (l,-)) may be unbounded and N{dt,dz) is its 
compensated Poisson random measure with N{dt 1 dz) := N(dt,dz) — v(dz)dt. The parameters b,a,j satisfy 
appropriate conditions (see Section [J) to ensure the well-definedness of this SDE. 

If an admissible control policy V = (ti, £i; T2, £2; • • •) is adopted, then X(t) evolves as 

O' 
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dX{t) = fj,(X(t-))dt + a(X(t-))dW(t) + [ j{X(t~), z)N(dt, dz) + V S(t - t;)&, 



(1.2) 



where S(-) denotes the Dirac delta function. With this given control, the associated total expected cost 
(objective function) is 



(r°° 00 \ 

J o e- rt /(A-(t))dt + 5^e-"'5(fi)J 



(1.3) 



> 
On 

The aim is to minimize the total cost over all admissible control policies, with the value function 

m 

(Np u(x) = MMV}. (1.4) 

&\ , HJB and Regularity. A heuristic derivation from the Dynamic Programming Principle shows that the 

value function (|1.4p is associated with the following Quasi- Variational-Inequality, or HJB, by 

max(£u -f,u- Mu) = in R n . (HJB) 
!-h ' where A4(p(x) is the so called minimal operator such that 

Mcp(x)= inf (cp(x + + B(0), (1-5) 
and C(p(x) is the partial integro-diffcrential operator 
&p(x) = — tr [A ■ D 2 <f(x)] — /i(x) ■ D(p(x) + np(x) + / [tp(x + j(x, z)) — (p(x) — j(x, z)Dip(x)] v(dz), (1.6) 

where the matrix A is given by A = {aij) nxn = \<j(x)a{x) T . Most recently [H] proves rigorously that indeed 
the value function is a continuous solution to (jHJBp in a viscosity sense. 
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Nevertheless, an important question remains: under what conditions is the value function a solution to 
the (|HJBp in a classical sense? Or, what is the degree of the smoothness (i.e. regularity property) for the 
value function in general? This is the focus of our paper. 

Regularity property has been one of the central topics in PDEs theory [31] [TBI [33]. Besides its obvious and 
natural theoretical interest, regularity study provides useful insight for numerical approximation. Controls of 
the impulse type by allowing discrete state space and fixed cost proves most desirable for application purpose. 
See PJJl [H C2J EH] for risk management, gHl[37] for real options, [1 [Ml EEE H HI HO] for transaction cost 
in portfolio management, [26[ [7J for insurance models, [35[ [6] for liquidity risk, and [271 ESI [8] for optimal 
control of exchange rates. Meanwhile, jump diffusions such as Levy processes, have been very popular in 
financial modeling. See for example [S EH [HI El [UJ [50l [3J . 

Combined, there is a growing interest and need to analyze impulse controls on jump diffusions. Unfor- 
tunately, impulse control is among the hardest to analyze and the regularity study for the associated HJB or 
the value function is largely open, except for some special and degenerate cases such as singular control and 
optimal stopping problems, see [351 33J [JJ51 [2] . One of the difficulty in establishing the regularity property lie 
in the non- linear, non-local operator M.u in Eqn (|HJB[) . Another difficulty is the partial integro-differential 
operator £ip(x) associated with the jump processes. For the special case when the controlled diffusion is 
without jumps, [3J established the regularity property by assuming that the control is bounded and non- 
negative with additional smoothness in the cost structure. Recently, [5D] applies the tricks of translating the 
regularity of the minimal operator in the action region into that of the PDEs in the continuation region. 
However, all these technique fail for a controlled jump diffusion component. The major issue is the addi- 
tional partial integro-differential operator. Moreover, (possibly infinite) jumps through the boundary might 
potentially reduces the degree of smoothness for the value function. 

Our work. This paper investigates the regularity of the value functions for the jump-diffusion models 
with impulse control. Building on the existence of the viscosity solution to the HJB for the value function 
[33] and the trick of [5D] for the non-local minimal operator, we focus on the partial integro-differential 
operator in the continuation region. There are two distinct cases: when the jump is driven by compound 
Poisson process, or equivalently when the Levy measure is finite, the analysis is fairly straightforward by 
the standard Schauder's estimate from PDEs, as in [20] . For the most interesting case of infinite Levy 
measure, the key is to combine the classical L p theory with the "bootstrap" argument to obtain regularity 
of the partial integro-differential operator. Finally, to deal with the regularity along the "free boundary", 
appropriate penalty function is devised. Surprisingly, despite the added possibly infinite jumps, we have here 
the same regularity as in the diffusion case, at least when the jump satisfies certain integrability conditions. 

2. Assumptions and Notations. We first specify the exact mathematical framework for our problem. 
Given a filtered and complete probability space (^Ft)t>o^) satisfying the usual conditions, we have a 

controlled jump diffusion process X(t) as defined above at (|1.2| . An admissible impulse control V consists of 
a sequence of stopping times t\,t%, . . . with respect to J- t and a corresponding sequence of R"-valued random 
variables £i,£2, ■ ■ ■ satisfying the conditions 

< n < r 2 < • • • < n < . . . , 
oo a.s. as i— > oo, 

6eJv 4) V»>1. 

The associated total expected cost (objective function) is given by (|1.4[) where / is the "running cost" , B is 
the "transaction cost" and r > is the discount factor. We assume that all the randomness comes from W 
and A, so that the filtration F = (J- t )t>o is generated by W and N. 

We next specify detailed conditions on the coefficients to ensure the existence and uniqueness of 
as well as the conditions on / and B in Ej2j 

Throughout this paper, we shall impose the following standing assumptions: 
(Al) Lipschitz conditions on (i : R" -> R",cr : R" -> R nxm ,j : M" xl' ^ R": there exist constants 
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Cft,C a > and a positive function Cj(-) G L 1 n L 2 (R', z/) such that 

- fi(y)\ < C^x ~ y\, 
\\a(x) - a(y)\\ < C„\x - y\, Vx, y G K", z G R z . (2.1) 

\j{x,z) - j(y, z)\ < Cj(z)\x-y\, 

Assume also that 

j(x, •) G L 1 ^; u) for every x G R". (2.2) 

(A2) Ellipticity: There exists a constant A > such that 

aij{x)ZiZi >A|£| 2 , V^CeM", (2.3) 

where the matrix A = (oyOnxn = ^ (J ( a; ) <7 ( a; ) T - 
(A3) Lipschitz condition on the running cost / > 0: there exists a constant Cf > such that 

l/OcJ-ZGOI^C/ls-tfl, Vx,yGlR". (2.4) 
(A4) Conditions on the transaction cost function B : M" — > R: 

inf B(£) = A" > 0, 

B G C(R"\{0}), (2 5) 

|B(£)|— > oo, as £|— > oo, and 
[B(a) + B(6) > 5(6 + 6) + A, 6 G R". 

(A5) r>2C^ + C 2 a + / RI C](z)v{dz). 

Assumption (Al) ensures the existence and uniqueness of solutions to (|1.1|) (Cf. Theorem 9.1, Chapter 
VI, [13] )■ The condition (|2.2[) seems essential to our approach, in particular to establish the continuity 
property of the operator / in Lemma [3~2l Readers are referred to [17] or [41] for a more detailed discussion 
of Levy processes and jump diffusions. 

In view of Assumption (Al) the following definitions for operators L,I make sense. 

Ltp(x) = - tr [A ■ D 2 (p(x)] - p(x) ■ Dtp(x) + r<p(x), (2.6) 
where p, = fj, — J Hi j(x, z)v{dz) and, for Lipschitz continuous functions ip, 



I(p(x)= / [(p(x+j{x,z))-(p(x)]v(dz) (2.7) 

We also adopt the following standard notations for function spaces: 

UC(M™) = space of all uniformly continuous functions on R ra , 
W k ' p (U) = space of all LP functions with /3-th weak partial 

derivatives belonging to L P ,V|/3| < k, 
Wq' p (U) = the closure, in W k ' p -norm, of smooth functions with compact support in U, 
W££(U) = {/ G W k ' p (U'),V compact U' C U}, 

f G C k {D) : sup { ^^^l < oo,y\p\ < k 

x,v£D ( \x y\ J 
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3. Preliminary Results. We first establish some preliminary results under the assumptions (Al)- 
(A5). 

Lemma 3.1. The value function u(-) defined by (|1.4[) is Lipschitz. 

Proof. Given an admissible control V, and two initial states X\,X2, denote by X l {t) the solution of 
(jl.ip . Apply ltd formula (for jump diffusions) (Theorem 5.1. Chapter II, to Y{t) = \Z(t)\ 2 , where 

Z(t) = X\t) - X 2 (t), 

dY(t) = 2Z(t) ■ [(^(A 1 ^)) - fi(X 2 (t)))dt + (criX 1 ^)) - o-(X 2 (t)))dW] 
+ {o{X\t)) o{X 2 {t)))(a{X\t)) - a(X 2 (t))) T dt 

- f 2(j(X 1 (t),z)-j(X 2 (t),z))Z(t-)is(dz)dt 

J-Rl 

+ f [(Z(t-)+j(X 1 (t),z)-j(X 2 (t),z)) 2 -\Z(t-)\ 2 }N(dt,dz) 

Integrating from to t, taking the expectation and then using Assumption (Al), we obtain 

EY(t) - (an - x 2 ) 2 < + Cl+ J C 2 {z) V {dz)^j j EY(s)ds, 

which implies that E|X x (f) - X 2 (t)\ < e ct \xi - x 2 \ by Gronwall's inequality, where C = 2(7^ + C 2 + 
J Ri C 2 (z)v{dz). Hence J Xl [V] — J X2 [V] < C u \x\ — x 2 \ by Assumptions (A3) and (A5), where 



r-[2C, + C 2 +J Ml C 2 (zMdz)] ■ 

By the arbitrariness of V, 

u(xi) < J X1 [V) < J X2 [V] + C u \xi - x 2 \ => u{x{) < u{x 2 ) + C u \xi - x 2 \. 

Exchanging the roles of xi,x 2 we get the desired result. □ 
Lemma 3.2. Iu e C(R"). 

Proof. Given x £ R", u(y + j(y, z)) — u(y) — > u(x + j(x, z)) — u(x), as y — > x, for any zEl'. Observe 
that if \y — x\ < 1, 

Hy + j(y,z))-u{y)\ <C u \j{y,z)\ <C«(\j(x,z)\ + Cj(z)\x-y\) < C u (\j{x,z)\ + Cj(z)), 

where C u is the Lipschitz constant of u. Since j(x,-) and Cj(-) are both i/-integrable, the dominated 
convergence theorem yields the desired result. □ 

For reference, we recall here Lemmas I3.3ll3.5l which were proved in |20j . 

Lemma 3.3 (Properties of M). 

1. M is concave: for any (pi,(p 2 G C(R n ) and < s < 1, 

M(sipi + (1 - s)ip 2 ) > sMipi + (1 - s)Mip 2 . 

2. M is increasing: for any ifi < (p 2 everywhere, 

Mfi < M<p 2 . 

3. M maps UC(W n ) into UC'(W l ) and maps a Lipschitz function to a Lipschitz function. In particular, 
Mu(-) is Lipschitz continuous. 

Lemma 3.4. u and Aiu defined as above satisfy u(x) < A4u(x) for all x £ R n . 
We define the continuation region C and the action region A as follows, 

C:={xe R" : u(x) < Mu(x)}, (3.1) 
A:={xe R" : u(x) = Mu(x)}. (3.2) 
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Then, C is open, and we have 

Lemma 3.5. Suppose x £ A, then 

(i) The set 

S(x) := {£ £ R" : Mu{x) = u(x + £) + B{£_)} 

is nonempty, i.e., the infimum is in fact a minimum. 

(ii) Moreover, for any £(x) £ S(x), we have 

u(x + £(x)) < Mu(x + £(x)) - K, 

in particular, 

x + ((x) e C. 



4. Viscosity Solutions. There arc different ways to define viscosity solutions. Let us begin with the 
most common one. 

Definition 4.1. A function u(-) £ UC(W n ) is called a viscosity subsolution (super solution, resp.) of 
(|HJB[) if whenever ip £ C 2 (M. n ),u — tp has a global maximum (minimum, resp.) at xq and u(xq) = <p(xo), 
we have 

msx{£(p(x ) - f(x ), <p(x ) - Mip(x a )} < (> resp.); (4.1) 

and u is called a viscosity solution of (jHJBI) if it is both a subsolution and a supersolution. 

Besides this standard definition of viscosity solutions, there are at least another two different (but 
equivalent) ones. The second way is to use semijets in stead of test functions. See, for instance, [T3] and [4*4"] . 
for more details. For the purpose of proving our regularity results in Section [5J we give a third definition 
below. The idea is that we impose "local" conditions (rather than global conditions as in Definition 14. ip on 
the test functions, and in the equation we only replace u by the test function ip in the "local" terms while 
still keep u in the "nonlocal" terms. The same definition (in different notation) and the proof of equivalence 
can be found in [35]. See also [3S][T] and [50] for a similar treatment. 

Definition 4.2. A function u{-) £ UC(W l ) is called a viscosity subsolution (supersolution, resp.) of 
(|HJB[) if whenever tp £ C 2 (R n ), u — ip has a local maximum (minimum, resp.) at xq and u(xq) = <p(xo), we 
have 

ma.x{Lip(xo) — f(xo) — Iu(xo),u(xq) — Mu(x n )} < (> resp.). (4.2) 

u is called a viscosity solution of (|HJB[) if it is both a subsolution and a supersolution. 
Theorem 4.3. The above two definitions of viscosity solutions are equivalent. 
Proof. See [35] Proposition 5.4]. □ 
We now have the following basic result. 

Theorem 4.4 ([44] [44]). The value function u(-) defined by (|1.4|) is a viscosity solution of (|HJB[) . 

This theorem was proved in [3T] Theorem 9.8] as well as [331 Theorem 4.20 in the sense of our Definition 
14.11 But when wc prove the regularity result below, we found it more convenient to use Definition 14.21 More 
precisely, by Theorem [331 and Theorem 14. 41 we can "identify" our value function u(-) with that of an impulse 
control problem of diffusion processes without jumps. 

Corollary 4.5. The value function u(-) is a viscosity solution of 

man{Lu(x) - f(x),u(x) - Mu(x)} = in R n , (4.3) 

where f{x) = f(x) + Iu(x). 

In 1441 the result was proved using an in principle smaller class of controls, the so-called 'Markov controls'. However, this 
restriction is unnecessary, as can be seen from the proofs of the analogous results in [48] or [25] . or |51j . 
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5. Regularity of Value Function. In this section we study the smoothness of the value function u, 
starting with the special case of a finite Levy measure. 

5.1. Special Case: u(M. ) < oo. Let us first consider the special case in which the Levy measure is 
finite, or equivalently, the jump diffusion X(-) is driven by a compound Poisson process. Then the operator 
/ enjoys the following nice property. 

Lemma 5.1. Suppose is(M. 1 ) < oo, then the operator I maps a Lipschitz function to a Lipschitz function. 

Proof. Suppose f{x) is Lipschitz with \(p(x) — if(y)\ < C v \x — y\ for any x,y £ R™, then 

\I(p(x) - I(p(y)\ < / \ip(x + j(x,z))-(p(y + j(y,z))\v(dz)+ / \tp(x) - <p(y)\v(dz) 
<C V ( [2\x -y\ + \j(x, z) - j(y, z)\] v{dz) 
< C v (2v(R l ) + Cj(z)u(dz)j \x - y\. 

So I(p is Lipschitz. □ 

Corollary 14.51 and Lemma 15.11 together imply the regularity of u in the continuation region. 

Lemma 5.2 (C 2 ' a -Regularity in C). Assume that a £ C 1 (IR") and v(R l ) < oo, then for any compact 
set D C C, the value function u(-) is in the Holder space C 2,a (D) for any a £ (0,1), and it is a classical 
solution of 

Cu — f(x) = in C. 



Proof. Note that u is a viscosity solution of (|4.3[) by Corollary 14. 5i and hence a viscosity solution of 
Lu — / = in C. On the other hand, / £ C a for any a < 1 by Lemma \5. II Classical Schauder estimates 
imply the desired results. (See the proof of Lemma 15.41 below for a similar argument.) □ 

Finally, an argument as in [20( §4] applies and yields the following 

Theorem 5.3 (W£' c p -Regularity). Assume that u(W) < oo and 

a £ C X > X {D) for any compact set D C E™. (5.1) 

Then for any bounded open set O C K™ and p < oo, we have u € W 2 p (0). 

5.2. More General Case: j(x, ■) S L 1 ^). Next, wc would like to remove the strong assumption that 
the Levy measure v is finite, and assume only our standing assumptions (A1)-(A5). Again, we first consider 
the regularity of u in the continuation region C, in which the linear elliptic PDE is satisfied. The difficulty 
is that we do not know lu is Lipschitz or even Holder continuous, but only continuous, by Lemma 13.21 
We cannot apply Schauder estimates at this stage, but the LP estimates give the following 
Lemma 5.4 (Wj^f -Regularity in C). Assume that a £ C 1 (M n ), then for any compact set D C C, the 
value function u(-) is in the Sobolev space W 2,P (D) for any p < oo, and it is a strong solutioi^ of 

Cu — f(x) = in C. 



Proof. Denote by / = / + lu, which is continuous by Lemma l3.2l Consider in any open ball B C C the 
following Dirichlct problem 

{ LW = l mB 1R (5.2) 
\w = u, on oB. 



2 A strong solution is a twice weakly differentiable function in the bounded domain that satisfies the equation almost 
everywhere. 
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Classical LP theory (Cf. [TSJ Corollary 9.18]) asserts that the Dirichlct problem (|5.2[) has a unique strong 
solution w G Wfo£(B) fl C(B) for any p < oo, since / £ C(B) and the boundary data u 6 C'(dB). Because 
a G C 1 (B),jl G C 0,1 (B) and / G C(B), this solution u; is in fact also a viscosity solution of (|5.2[) by [Ml 
Theorem 2]. 

On the other hand, it is also a viscosity solution of (|5.2[) by Corollary 14.51 Therefore, iu = u in B 
by classical uniqueness results of viscosity solutions to a linear elliptic PDE in a bounded domain(Cf. [T31 
Theorem 3.3]). Hence u G W^{B) n C{B). 

Finally any compact set D C C can be covered by finitely many balls {B r (a?fe)}^L 1 of radius r < 
i dist(.D, <9C). Let -B = B2 r (xk) C C in the above argument, then u is in W 2 ' p (B r (xk)) for all k and also in 
W 2 ' P (D). 

With more regularity of u in the continuation region C, we can use the "bootstrap argument" to obtain 
further regularity of Iu (and hence it) in C. 

Theorem 5.5 (C 2 '"-Regularity in C). Assume that a G C 1 (R"), then for any corn-pact set D C C, the 
value function u(-) is in the Holder space C 2,a (D) for any a G (0, 1), and it is a classical solution of 

Cu — fix) = in C. 



Proof. The key step in the proof is to show Iu G C a (D) for any compact D C C. 

Take a compact set D' such that D C D' C C and 5 := dist(D, <9L>') < 1. Then by Lemma HT2l 
u G W 2 > V {D') for any p < oo. By Sobolev imbedding, u G C l a (D') for all a G (0, 1). Define the set 

E :={z eR l : \j(x,z)\ < 6,Vx G £>}. 

Then for z e E c =S}\E, there is x G D such that 

5< \j(x,z)\ < \j(0,z)\+C 3 (z)\x\ < \j(0,z)\ +000,(2), 

where Co = max{|x| : x G D} is a constant. So |j(0, z)| > <5/2 or Cj(z) > S/(2Cd) and 

K£ c ) < 7 / U(o^)IK^) + ^ / Cj{zMdz) < 00. 

For any xi, X2 G D, 

\Iu(x x ) - Iu(x 2 )\ < [ \[u(x 1 + j{x 1 , z)) - u(xi)] - [u(x 2 + j(x 2 , z)) - u(x 2 )]v(dz) 

\u(xi +i(xi, 2)) - u(x 2 +j(x 2 ,z))\ + \u(xi) - u(x 2 )\v(dz) 



E 



IE" 

< \Du(x! + sj(x!,z)) ■ j(xi,z) - Du(x 2 + sj(x 2 ,z)) ■ j(x 2 ,z)\ds v(dz) 

J E JO 

+ \xi-x 2 \ C u {2 + C j {z))v(dz). 

J E c 

< \Du(xi + sj(xi, z)) — Du(x 2 + sj(x 2 , z))\ ■ \j{x\, z)\ds v(dz) 

J E JO 

\Du(x 2 + sj(x 2 ,z))\ ■ \j(x 1 ,z) -j(x 2 ,z)\ds vidz) 

IE Jo 

+ |xi-x 2 | / C u {2 + C j (z))u(dz). (5.3) 

JE" 

Note that xi + sj(xi,z), x 2 + sj(x 2 , z) G D' for all < s < l,z G E and that Du G C a (D'). Thus the first 
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integral in (|5.3p can be estimated by 

\\ Du \\c«(D>) J (^J \xi- x 2 + s(j(xx,z) - j{x 2 ,z))\ a ds S j\j{x 1 ,z)\v{dz) 
<\\Du\\c~(d>) I \x 1 -x 2 \ a (l + C j (z) a )\j(x 1 ,z)Hdz) 



/ (1 + C j (z) a )\j(x,z)\l {z .\ j ( x , z)l< i } u(dz) 

\j(x,z)\u(dz)+ ( Cj(z) 2 v(dz)] ( \j{x,z)\^l {z: \ j(XtZ)l<1} v(dz) 



< C\\x\ - x 2 \ a , 

for some constant C\ > independent of x±,x 2 , because by Holder's inequality. 



< 



< 

< / \j(x,z)\v(dz)+[ I C 3 {z) 2 v(dz)\ ( / \j(x,z)\v(dz) 

which is a continuous function in x and has a maximum on D independent of X\, x 2 . 

The second term in (|5.3p can be majored by ||Du||j,oom) J Ri Cj(z)v{dz)\x\ — x 2 \ =: C 2 \x\ — x 2 \ and the 
third term is majored by C u [2i>(E c ) + J Ri Cj(z)v(dz)\\x\ — x 2 \ =: C3|a;i — x 2 \. Thus, putting all three terms 
in (|5.3p together, 

|/(flfi) - I(x 2 )\ < d\ Xl - x 2 \ a + (C 2 + C 3 )|x! - x 2 \ < C\xi - x 2 \ a , 

where Ci,C 2 , C3 and C = C\ + (C 2 + C3)(diam D) l ~ a are constants independent of xi, x 2 G D. This proves 
that Iu G C a (D). 

Finally, we can repeat a similar argument in the proof of Lcmma l5.2l This time we know / G C a (D), thus 
the solution of (JSHJ) is in fact in C 2 a {D) by Schauder estimates (PH Theorem 6.13]). Thus, u G C 2 a {D) 
for compact D C C, and u is a classical solution of Cu — / = in C. □ 

Parallel to [20], once we have C 2 ' a regularity, we are able to obtain W 2,P (D) regularity for any compact 
set Dei". 

Theorem 5.6 (W£' c p -Regularity). A ssume that 

a G C ls (D) for any compact set D C R™. (5.4) 

Then for any bounded open set O C K™ and p < 00, we have u G W 2,p {0). In particular, u G C 1 (R n ) by 
Sobolev embedding. 

The proof is given in the next section. It follows similar lines to that of [201 Theorem 4.2], involving 
regularity of an associated optimal stopping problem. 

6. Proof of Theorem 15.61 To study the regularity of the value function for the impulse control 
problem, we need to investigate the related optimal stopping problem. More precisely, we shall obtain the 
regularity of solutions for the HJB equation associated with this optimal stopping problem. 

Theorem 6.1. Suppose O is a bounded open set in W 1 with smooth boundary. Assume 

CHjGC^iO), meC^iO), r>0, feC(0), (6.1) 

C > c|£| 2 ,Vx,£ G O, for some c> 0, (6.2) 

g G C(0), g > on dO. (6.3) 
Assume also that there exist a sequence of functions {g 6 } e >0 and a constant M > satisfying 

IV G C 2 (0) n C(0), _Lg" > —M m O, 
I g e — > g uniformly in O. 



If v G C{0) is a viscosity solution of 

fmax{Lv -/>-#} = in O, 

\v = on dO, 1 ' ' 

then v S W 2 ' p (0). 

Remark 1. Note that Assumption (|6.4p is trivially satisfied if g £ C 2 (0). However, later we will apply 
this theorem to g = Mu, which is not necessarily in C 2 (0). In applications, g e can be taken as the usual 
mollification of g, or its slight modification (which may only be in C (O) fl C(0) but not in C 2 (0), as in 
Corollaru \6.S\ below). 

As a corollary of Theorem 16.11 we obtain local W 2,p (n < p < oo) regularity of continuous viscosity 
solutions of 

m&x{Lv - f,v -g] = in R". (6.6) 

Corollary 6.2. Assume that f G C(R n ), a i3 G C^(R n ), and (i, a and g are Lipschitz in M. n . Assume 
also that for any bounded open set O C M. n with smooth boundary, there are constants ( maybe depending on 
O) c > and M such that (|6.2| and (|6.4|) are satisfied. 

If v G C(R") is a viscosity solution of (|6.6[) . then v G W 2,p (0) for any O C R n with smooth boundary 
and any 1 < p < oo, and hence also in C 1 . 

We defer the proofs of Theorem 16.11 and its corollary to the appendix, and focus now on proving our 
main theorem using the above corollary. 

Proof. [Proof of Theorem 15.6] Given any bounded open set O with smooth boundary, we denote by C 
(A 1 , resp.) the restriction of the continuation (action, resp.) region within O. Then there exists an open 
ball O' D O such that for any x G O, u(x + £) + B(£) < Aiu{x) + 1 implies x + £ G O' . Because in this case, 

B{0 < Mu{x) -u(x + {;) + 1 < Mu{x) + 1 < supMa + 1 < oo. 

o 

But £?(£)— > oo as » oo, which implies that all such £ are bounded uniformly. 
Now we define the set 

D;= LeO':u(y)<Mu(y)-^Y (6.7) 

Clearly, D is compact and D C C. From Lemma [521 

u G C 2 - a (D). 

For any x € O, take a minimizing sequence such that u(x + £,k)+B(£,k) — * Mu{x). Then C O'. 
Extract a convergent subsequence (still denoted by {£fc}) converging to £*. Because -B(£) + £?(£') > K + 

Mu(i) = inf {w(a; + £ fe + r?) + 5(6 + ry)} 

77GR" 

< inf {u(x + ^ +T1 ) + B(n)} + - # 

= X M (x + 6) + S(a)-/f 

= JHu(i + 4) - u(a: + &) + [«(a; + &) + 5(4)] - 
Passing to the limit k— > 00, we obtain 

u(x + £*)- Mu(x + C) < ~K. 

In particular, y := x + £* G -D. 



On the other hand, since u — A4u is uniformly continuous on O' , there exists po > such that 
\y - y'\ <Po^ Hy') - Mu{y') - (u(y) - Mu{y))\ < -. 
Hence, for all p G (0, po],^ G [—1, 1] and unit vector \ £ 

y = x + £* G D, y' = y + X PX G D, 

because u(y') - Mu(y') < u(y) - Mu(y) + -f < --§■• 

Since Mu(x ± px) < u(x ± px + £fc) + B(£k) for all fc, 

Mu(x + px) + Mu(x - px) - 2A4u(x) 
<u(x- + px + &) + u(a: - + &) + 2.B(&) - 2A4u(*) 
^u(y + px)+u(y-px)-2w(?/), fc-> oo, 



and hence the second order difference quotient at x 



-^r[Mu(x + px) + Mu(x - px) - 2Mu(x)] 
P 

<~2 I M (Z/ + PX) + u{y - PX) - 2u(y)] 
P 

1 f 1 

=r7 / l( Du (y + x px) - Du (y - A/?x) • x^a 

\P\ Jo 

<c D , 

where Cd = sup \D 2 u(x)\ < \\u\\ C2 , a my 

For simplicity, denote by g = A4u and g E its mollification. For any xq G O, suppose Bg(xo) C 0, then 
for any e G (0, §), p G (0, p A f ), and a unit vector x G R", 

-j [fl^o + /°%) + 9 £ {xo - PX) ~ 2ff £ (-To)] 
P 

= ~o / [5(^0 - z + px) + g{x -z-px)~ 2g(x - z)]rf{z) dz 

<C D f V 6 {z)dz = C D . 
Jb e (o) 

Sending p — > we get 

X T #y(a;o)x < C D . 

Hence, 

tr(a(x )a(x ) T D 2 g £ (xo)) = tr(a T (x )D 2 g £ (x )cr(x )) 



fc 



<C, 

where Cfc is the fc-th column of the matrix a, <Tij is the (i, j)-th element of cr, and the last inequality is due 
to continuity of a. 
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Note that |g E (cco)| + \Dg e (xo)\ < \\g\\wi-°°(0) an d n(x) bounded, we deduce 

Lg E (x ) = -^ti(a(x )a(x a ) T D 2 g e (x a )) - fx(x ) ■ Dg E (x a ) + rg E (x a ) > —M, 

where the constant M is independent of xq. 

Finally, recall that u is a viscosity of (|4.3[) by Corollary 14. 51 We can apply Corollary 16 . 2 1 with / replaced 
by / = / + Iu G CiW 1 ) and g = Mu and conclude that for any 1 < p < oo, 

u g T47 2 ' p (0). 

□ 
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Appendix A. Proof of Theorem l6.ll A standard technique to get regularity is to consider a sequence 
of penalized problems. For this, let [3 e denote a sequence of smooth functions satisfying 

{P £ (t) — ► oo, as e -> 0,t > 0; 
(3 £ {t) -> 0, as e -> 0,t < 0; 
0<&'(t) <Lo(e)-\Vt; 
/3 e (0) = 0,/3 e >-l. 

Here tu(-) is the modulus of continuity for the convergence g £ — > g, i.e., 

w{8) := sup ||fif e — 5llc7(o)- 

Thus, oj(e) — > as e — > 0. For a given e > 0, the graph of (3 E is shown in Figure [A~T1 
We approximate the Dirichlct problem (|6.5p with the following penalizing problems: 

Lv E + /3 E (v e - g e ) = f inO, 
v e = on dO 

By a standard fixed point argument, (|A.2[) has a unique solution in W 2,p (0) n W 1 ' p (O), for any 1 < p < oo. 
Moreover, we have the following estimates: 

Lemma A.l. Under the same assumptions as Theorem \6.1[ there exists a constant C independent of e, 
such that 



(A.l) 



(A.2) 



v e \\w 2 -p(o) < C. 
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Proof. The goal is to show that 

m(v £ - 9 £ )\\l~(o)<C. (A.3) 

Clearly, to get (|A.3j) . it suffices to show 

f3 e (v £ - g £ ) < C, (A.4) 

since j3 e > — 1 by construction. 

Consider the point xq at which the maximum of f3 e {v e — g e ) occurs. 
Case 1. If xq £ dO, then v s (xq) — < g(xo). Thus 

v £ (x ) - g e (x ) < g(x ) - g £ {x ) < w(e). 

Because < (3 £ ' < w(e)- 1 and (3 e (0) = 0, we have (3 e (v £ - g e ) < 0, if v £ - g £ < 0; or < /? e (V - g £ ) < 1 if 
v £ — g e > 0. In either case, we have (3 e (v e — g e ) < 1 at xq. 

Case 2. If xq £ O, then v e — g e also attains its maximum there, since (3 e ' > 0. By Bony's maximum 
principle (see [5] or [3"1]). 

ess lim sup{-aj,(w e - g 5 )^.} > 0. 

X — >Xq 

Thus, either v £ — g £ < at xq, which implies (3 e {v e — g e ) < 0, or, 

ess lim supLv 6 > Lg e {xo) > —A/, 

x — yxQ 

due to (|6.4p . By continuity, 

/? e (u e - 5 e )(x ) = f(x Q ) - ess lim supLt; 6 < f(x ) + M < C. 

X — >Xq 

Combing these two cases, we obtain (|A.4[) . Thus, we proved (|A.3[) . and from the PDE (|A.2[) we have 
||^^ e ||i°°(0) < C independently of e. 

Finally, the Calderon-Zygmund estimate implies the desired result. □ 
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Thanks to Lemma IA.11 we can extract a subsequence of v £ , denoted by v £k , such that 



I v £k ->> v weakly in W 2 > p (0), 
1 v £k — > u uniformly in O. 

Due to the stability result of viscosity solutions (see, for instance. [121 Theorem 8.3]), this limit function v is 
in fact a viscosity solution of (|6.5[) . 

Lemma A. 2. The limit function v is a viscosity solution of (|6.5[) . 

Proof. (Subsolution) Suppose <f)(x) is a smooth test function and v — (f> has a local maximum at xo E 
with w(xo) = 4>(x$ ). We want to show that 

max{L0 — /, cf> — g} < at xq. 

Without loss of generality, we may assume Xq is a strict local maximum because otherwise we can replace 
<f>(x) by 4>(x) — \x — Xo| 4 and prove the same result. In this case, for any open ball with radius S > and 
centered at xq, denoted by Bg(xo), v £k —cj> has a local maximum x k E Bs{xq) for /c sufficiently large, because 
v £k converges to v uniformly in B$(xq) and v(xq) — 4>( x o) > 1 ^XdB s (x ){^ — 4>)- Let S go to zero, and we 
obtain 

Xk is a local maximum of i> £fc — </>, and lim Xfc = xo- 

By the maximum principle, ess lim sup{— ciij(v £k — (j>) XiXj } > at Xk- However, v £k satisfies (|A.2I) . 
thus, 

L(j>{xk) < ess lim supLv £k — r(v £k (x k ) — 0(#fc)) 

= /(x fe ) - &> e * - g Ek )(x k ) - r(v £k (x k ) - <f>{x k )). (A.5) 
Note that v £k — g £k — > C — g and i> £fc — <j) — > v — <f> locally uniformly, we have 

L £k (x k ) -g £k (x k ) ->v(x ) -g(x ), ^ ^ ^ 
[v £k [x k j - cf>{x k ) -> v(x ) - <p(x ) = 0, 

Sending k— > oo in (|A.5[) . we obtain 

4> — g = v — g < at xq ■ 
(Otherwise, Lcf) < — oo at xq since /? e (t)— » oo as s — > if t > 0). Moreover, 

L(j) < / at Xq, 

since /3 £ (t) — ► as e — ► if t < 0. We have proved that v is a viscosity subsolution. 

(Supersolution) Similarly, suppose (j>(x) is a smooth test function and v — <fi has a local strict minimum 
at xq with u(xo) = 0(a;o). We want to show that 

max{L0 — /, cf> — g} > at xq. 

For the same reason as above, we can take a sequence {x k } so that 

x k is a local minimum of v £k — 4>, x k — > xo, 

and (|A.6[) still holds. And again by maximum principle, at x k , 

Lcf) > f - (3 ek (v 6k -g £k )-r(v £k - 0). (A.7) 
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If <P — 9 > at xo, then we have the desired inequality. Otherwise, 0(xo) — g{xo) = — 2i/ for some f > 0. 
So v Sh {xk) — g £k (xk) < ^ v < for fc sufficiently large by (|A.6[) . and hence sending oo in (|A.7[) yields 

L<j) > f at xq. 

Thus, v is a supersolution. 

Finally, the boundary condition is satisfied since v = on <90. □ 

Proof. [Proof of Theorem 16. 1| Because of the above Lemma fA. 21 and the uniqueness of viscosity solutions 
for pTlBj) . we conclude that v = v G M /2 ' P (C). □ 

Proof. [Proof of Corollary 16.21 ] To apply Theorem 16.11 we need to subtract from v a function that has 
the same boundary value on dO. Let w be the unique solution of the Dirichlct problem 

Lw = in O, 
w = v on dO. 

Then w £ C^"(0) D C(0). Thus, vq := v — w is a viscosity solution of 

max{Lu - /, v - g} = in 0, w = on dO, (A.8) 

where g = g-w. Then £ = .g- v > on <9(9. Take g £ = g e -w G C" 2 (O)nC"(0), satisfying L<f = L ff e > — M 
in 0. All the other conditions of Theorem [63] arc easily verified. So we have vq G W 2,p (0) and v = vq + w € 
W?*(0). But since O C W 1 is arbitrary, we have v G W 2 'P(£>) for any Ocf'.D 
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